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Abstract 

Midpoint subdivision generalizes the Lane-Riesenfeld algorithm for 
uniform tensor product splines and can also be applied to non regular 
meshes. For example, midpoint subdivision of degree 2 is a specific 
Doo-Sabin algorithm and midpoint subdivision of degree 3 is a specific 
Catmull-Clark algorithm. In 2001, Zorin and Schroder were able to 
prove C 1 -continuity for midpoint subdivision surfaces analytically up 
to degree 9. Here, we develop general analysis tools to show that the 
limiting surfaces under midpoint subdivision of any degree > 2 are 
C 1 -continuous at their extraordinary points. 
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1 Introduction 

The midpoint subdivision schemes form a class of subdivision schemes for 
arbitrary two-manifold meshes. The midpoint subdivision scheme of degree 
n G N is given by the operator 

M n = A n ~ l R , 

where R and A are the refinement and averaging operators, respectively 
Refining a mesh M by R means to connect the center of each face of M with 
the midpoints of all its edges, which results in the quadrilateral mesh RM. 
Averaging M means to connect the centers of all adjacent faces, which results 
in the mesh AM, see Figure (1.1). 




(1.1) Figure. Masks of R (left and middle) and A (middle) and an example 
for M 2 M (right). 

The midpoint subdivision schemes of degree 2, 3 and 4 are specific in- 
stances of the well-known Doo-Sabin, Catmull-Clark and Qu algorithms 
[DS78, CC78, Qu90]. 

The valence m of a vertex or a face is the number of incident or surround- 
ing edges. An inner vertex or face is called regular if m — 4 and irregular 
or extraordinary otherwise. After subdividing a mesh by M n , the mesh has 
no irregular vertices or faces if n is even or odd, respectively. Subdividing 
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the mesh further does not increase the number of extraordinary vertices and 
faces and extraordinary elements become more and more isolated. 

Therefore, and since midpoint subdivision generalizes the Lane-Riesenfeld 
algorithm, midpoint subdivision surfaces are spline surfaces except for finitely 
many extraordinary points. At these points, the smoothness analysis is 
complicated. Using Reif's (^-criterion [Rei95, Theorem 3.6] and interval 
arithmetic, Zorin and Schroder [ZS01] showed C^-smoothness for degrees 
ra = 2, ...,9. 

Naturally, this numeric approach is limited to a finite number of degrees. 
Here, we develop a geometric framework that allows us to prove (^-continuity 
of midpoint subdivision surfaces of all degrees. Moreover, we think that this 
framework may conceptually be useful for other classes of (simple) subdivi- 
sion schemes and that it provides - for the first time - an analysis toolbox 
for a complete class of subdivision schemes. 

In Section 2, we discuss the basic topological dependencies between the 
vertices of a mesh before and after a subdivision step. These dependencies 
correspond to a certain block structure of the subdivision matrix with a block 
possessing a strictly positive power. This particular block represents the sub- 
division operator restricted to a certain central part of the mesh, which we 
call a core mesh. In Section 3, definitions and properties of symmetric meshes 
and symmetric subdivision schemes are introduced. For symmetric meshes, 
in Section 4, we define a partial order based on particular coordinate sys- 
tems. In Section 5, we show that subdividing specific symmetric grid-like 
core meshes results, in the limit, under normalization in a symmetric mesh 
whose regular vertices do not coincide with the center and show that it is 
an eigenvector of the subdivision matrix restricted to core meshes. In Sec- 
tion 6, we use the partial ordering introduced in Section 4 to compare the 
eigenvector constructed in Section 5 to any other eigenvector and show its 
subdominance. Hence, in Sections 5 and 6, we restrict the analysis to core 
meshes for which the subdivision matrix has a strictly positive power. Then, 
in Section 7, we extend the analysis to larger meshes using the particular 
block structure of the subdivision matrix. In Section 8, we show that sub- 
dominant eigenvectors define a regular surjective characteristic map, which 
concludes the proof that generic midpoint subdivision surfaces of any degree 
n(> 2) are (^-continuous. This result is stated in Theorem (8.4) and the 
entire paper consists of its proof which is composed of 19 lemmata, theorems 
and corollaries whose ramified interdependencies are depicted in Figure (1.2). 
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L: Lemma 
T: Theorem 
C: Corollary > 



(2.2) L 



Section 2 
C°-Property 




Section 3 



Section 7 



Section 8 
C 1 -Property 



(1.2) Figure. Interdependencies of lemmata, theorems and corollaries in 
this paper. 



2 Rings and ringnets 

For the smoothness analysis of midpoint subdivision surfaces at extraordinary 
points, it is sufficient for odd n to consider a mesh with only one irregular 
vertex and for even n to consider a mesh with only one irregular face. We 
will assume this throughout the paper. These simple meshes are illustrated 
in Figure (2.1) and are called ringnets. 

The 0-ring of a ringnet is formed by the irregular vertex or the irregular 
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n U TZ X 1Z U Tlx U ft 2 

□ e ft o e {^} = n nefix o e n 2 

(2.1) Figure. Example of rings and ringnets: a 1-ringnet with an irregular 
face of valence 5 (left) and a 1-ringnet with an extraordinary vertex of valence 
5 (right). 



face, the 1-ring by the vertices of the adjoining faces that do not belong to 
the 0-ring. The 2-ring, 3-ring, etc. are defined similarly. The rings 0,1,...,/ 
constitute an l-ringnet or an l-net for short. A ringnet Af is understood as 
a vector whose coordinates represent the vertices of the ringnet or as the 
set of its vertices, depending on the context. Following other authors, we 
call a ringnet with an extraordinary vertex primal and a ringnet with an 
extraordinary face dual. 

In this paper, n always denotes the degree of the midpoint subdivision 
scheme and m the valence of a vertex or of a face. In the following, let Af be 
a sufficiently large ringnet and let Af^ be the subdivided ringnet M^Af. If 
Af is primal we require n to be odd and otherwise to be even. The i-th ring 
is denoted by Afi, and the subnet built from Afi,Afi+i, . . . ,Afj is denoted by 
Afi..j. In particular, Afo..i is called an i-ringnet or i-net or i-mesh for short. 

In the rest of this section, we study on which vertices of Af the vertices 
of Af' : = Af^ depend and derive the minimal size of a ringnet required to 
examine the eigenvalues of the subdivision matrix. 

(2.2) Lemma. (The influence of vertices under subdivision) 

We say that Afi influences Afj and denote this by Afi ~* Afj if every vertex in 
Afi influences some vertex in Afj and if additionally all vertices in Afj depend 
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on Mi. This is the case if and only if 



2% 



n + 1 



< j < 2i + 



n + 1 



Consequently, Mq j is completely determined byMo..j, i. e., 

Mi 3 = {M n Mo.. j V.j 

for all j > oj := [ rk Y~\ ■ Furthermore, for i < u, it follows that 
M {0) ^ M {1 \ % M {2 K M {3 K 5& • • • ^ A r(r i °s 2 (-+Dl) 

and since Afo or (RAf) consists of only one vertex, every vertex in Mi influ- 
ences every vertex in Nq%, where k > [log 2 (cj + 1)] . 

Proof. The lemma can be shown via induction if one observes first that 

Mi £ (A 2 M)j if and only if i - 1 < j < i + 1 , 
secondly for a primal M that 

Mi 4 (BM) j if and only if 2i - 1 < j < 2i + 1 , 
and finally for a dual M that 

Mi (ARM)j if and only if 2i - 1 < j < 2i + 2 . 

□ 

The w-mesh An w with 



w = 



n 



is called the core mesh oiM (with respect to M n ). According to Lemma (2.2), 
the core mesh of M is the largest sub-ringnet of M with the property that 
each of its vertices influences every vertex in the core mesh after several 
iterations of subdivision. 

The meshes M {i) = M l n M converge to a piecewise polynomial surface. 
Each regular subnet of (n+ 1) x (n + 1) vertices of M^ defines a polynomial 
patch of this surface and all these patches form a spline surface ring Sj with 
an m-sided hole. The difference surfaces Sj\sj_i are smaller spline rings 
consisting of 3m [|J 2 polynomial patches determined by Mq} p with 



(2.3) p = 



-n 



see [Pra98]. The operator M n restricted to p-nets is represented by a quadratic 
matrix S p , called the subdivision matrix. 
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(2.4) Theorem. (C°-property of M n ) 

The subdivision surfaces generated by M n are C° for all n > 1 . 

Proof. Since the subdivision matrix S p is stochastic, 1 is the dominant eigen- 
value of S p . Using Lemma (2.2), we obtain that every vertex in J\f influences 
all vertices in S p Afo.. p . This implies that S 2 p has a positive column. According 
to [MP89, Theorem 2.1], the sequence (S l p c) converges to a multiple of the 
vector [1 ... I] 1 as % — > oo for all real vectors c. Therefore, the only dominant 
eigenvalue of S p is 1 and has algebraic multiplicity 1. 

Hence, the difference surfaces Sj\sj_i converge to a point, see [Rei95, 
Theorem 3.2] and [Che05, Remark 35 on Page 35], from which it follows that 
the surfaces generated by M n are continuous. □ 



3 Symmetric ringnets 

In this section, we introduce grid meshes and show that reflection symmetric 
eigennets of M n have real eigenvalues. 

(3.1) Definition. (Grid mesh) 

A primal grid mesh of valence m and frequency f is a planar primal ringnet 
with the complex vertices 

g\. = ie^if/m + Je ^(m)//m i,j>0, le Z m , % = V^l . 
A dual grid mesh of valence m and frequency f consists of the vertices 

h \j = + g{j_i + slu + slj), hi > i, z e z m , 

see Figure (3.2). For I fixed, the vertices g\- or h'- with (i,j) ^ (0,0) of 
a grid mesh J\f build the l-th segment of N '. The vertices g°- or h°- with 
i > j and ^ (0,0) constitute the first half segment of Af, denoted by 
H(Af). We call the vertices g l 0j = g^J 1 , g l u and h l u with i,j > 1 spoke 
vertices and call the vertices g^- or h' - with i ^ j and i,j>l inner vertices 
of the l-th segment. The segment angle of M is ip = 2irf/m. The half-line 
from the center g l 00 through g' 10 or through g l u is called the l-th spoke or the 
(I + \)-th spoke, denoted by Si{M) and S l+ i{H), or Si and S l+ i for short, 
see Figure (3.3). 

Topologically, any ringnet M. is equivalent to a grid mesh M . Therefore 
we use the same indices for equivalent vertices and denote the vertices of M. 
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(3.2) Figure. Primal grid mesh (left) and dual grid mesh (right). 



S, 



'S3 = So. 5 



f 




(3.3) Figure. The spokes of a ringnet with m = 5 and frequency f = 2 (left) 
and f — 4 (right). The segment angles are \ix (left) and §7r(= — fvrj (right). 



(3.4) Definition. (Symmetric ringnet) 

.A ringnet of valence m with the vertices p' ■ m C zs ca//ed rotation symmetric 
u^/i frequency f , if 

l J2nf/m _ l+l 
Vij^ Vij ■ 

A ringnet Af in the complex plane is called reflection symmetric if its permu- 



tation J\f consisting of the points p[j = p,-' 



(m-l)-l 



'31 



equals the conjugate ringnet 



Af, i. e., 



Af = Af . 

A rotation and reflection symmetric ringnet is called symmetric. 

It is easy to see that a symmetric ringnet is geometrically symmetric with 
respect to each spoke. Hence, we obtain 

(3.5) Lemma. (Eigenvalue of a reflection symmetric eigennet) 

If N is a reflection symmetric eigennet of M n with eigenvalue A and with 
j + 1 rings, j > oj, i. e., (M n Af)o..j — XAf , then A is real. 

Proof. Since 

XAT = (M n Af) ..j = (MS)o..j = (M n Af) ..j = A Af = \Af , 
we get A = A. □ 

(3.6) Lemma. (Positional invariance) 

Let Af be a symmetric ringnet with segment angle if G (0, n] . If the first half 
segment H(Af) of Af lies in the cone C{Af) spanned by the spokes S and S , 5 
ofAf, i. e., 

C(Af ) = convex hull of S (Af) and S , 5 (Af), 

then the subdivided ringnets RAf , AAf and consequently M n Af have the same 
property. By symmetry, the similar statement also holds for all half and, 
hence, for all full segments. 

Proof. The vertices of H(RJ\f) and H(AAf) are convex combinations of H(Af)U 
{0} or lie on So and S0.5 since R and A preserve symmetry. □ 

4 Comparing ringnets 

In this section, we partially order symmetric ringnets by comparing the co- 
ordinates of all vertices in the first half segment of any two meshes. Since 
we need that orthogonal projections into the spokes preserve this order, we 
choose the coordinate axes perpendicular to the spokes. This partial ordering 
is preserved by midpoint subdivision. Moreover, if Afi ^ A2 for two distinct 
core meshes, sufficiently many subdivision steps lead to the strict inequality 
M k Af\ > .U^A'j. where M is the operator M n restricted to core meshes, see 
Lemma (4.4). This will be crucial to prove that subdividing a grid mesh leads 
to a subdominant eigennet and in Lemma (4.6), we show that this eigennet 
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has no zero control points except for its extraordinary point if the net was 
primal. 

The coordinate system K of a planar symmetric ringnet M with segment 
angle (p G (0, 2n) has the basis vectors 



cos 9 
sin 6* 



and 




1 



, where 9 = (p/2 — n/2, 
which are orthogonal to the spokes 60.5 and So, see Figure (4.1). 

1 So. 5 



y 




Si 







x 



Sn 



(4.1) Figure. Examples ofKforcp/2 = Z(S , S . 5 ) G (0, f ) (left) and cp/2 
ASo, S0.5) = f (right). 



The midpoint of two vertices whose positions are symmetric with respect 
to So or S0.5 is the orthogonal projection of these vertices into So or S0.5, 
respectively. Thus, the two vertices and their midpoint have the same x or 
y coordinates in K, respectively. This property is crucial in our proofs, see 
Equation (4.8), Figures (4.5) and (4.9). 

A point p with coordinates pi,P2 in K is called positive if p\ > and 
P2 > 0, and it is called non-negative if p± > and p 2 > 0. We denote this by 
p > and p > 0, respectively. 

A symmetric ringnet M. is called positive or non-negative if all points of 
H(M) are positive or non-negative, respectively. This is denoted by M > 
or M > 0. 

(4.2) Definition. (Comparison of two symmetric ringnets) 

Let M\ and M2 be two symmetric ringnets with the same topological structure 
and coordinate system K. We call M\ > A2 if H(N\) > H(J\f 2 ) holds point- 
wise in K. The relations >, <, < and = between two symmetric ringnets can 
be defined similarly. 
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In order to characterize the size of a ringnet in K, we introduce a norm. 
(4.3) Definition. (Norm, MAX, MIN) 

Let M be a symmetric ringnet with coordinate system K. Let p be a vertex 
with coordinates pi and p 2 . Then we define 

||p|| := max{|pi|, \p 2 \} , 
MAX(AT) := \\Af\\ := max ||q|| 

qS-ff(Af) 

and 



MIN(jV) := min llqll . 
q eff(A0 

(4.4) Lemma. (Comparable properties between two symmetric ring- 
nets) 

Let M\ and J\f 2 be symmetric ringnets with segment angles ip^ e (0, n) and 
tp 2 € [y?i,7r], respectively, such that H{Ni) lies in C(A/i). //we rotate M\ by 
(ifi2 — ty?i)/4 degrees counterclockwise, we obtain the following statements with 
respect to the coordinate system 0/A/2 under the assumption Af± > A/2' 

r«; m„m > M n u 2 . 

(b) IfO < (A/i)o..u> 7^ (A/2)o..w> ^en ; /or any j > i/iere exzsfo an a G (0, 1) 

a(M£A/i) .j > (M*M)aj 

/or a// sufficiently large k. 

Proof. Because A/i is symmetric, H(M n Af i ) can be generated from H(J\fi) U 
{0} by calculating midpoints and by projecting points into the spokes So(A/i) 
and <So.5(A/i). Averaging preserves inequalities and the same holds for ortho- 
gonal projections into the spokes as we see for two corresponding vertices 

P i = [x i yi] t = i>k j (W, 1 = 1,2. 

Let be their projections into S (Afi) and let r^ be their projections into 
<So.5(A/i), see Figure (4.5). Since p x e C(A/i), we obtain 

Pi > P 2 =>• Qi > Q2 and n > r 2 , 
by which we derive Statement (a). 
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(4.5) Figure. Projection into the spokes. 

Moreover, since ipi < ir, we get 

< Pi > p 2 and X\ > x 2 =>- q x > q 2 , 
< Pi > p 2 and yi> y 2 ^ r 1 > r 2 , 

see Figure (4.5). We have assumed Af± > Af 2 and < p x > p 2 ^ p 1 for 
some points p i = p° 1 ; 2 ((A/i)o..u;), i = 1,2. According to Lemma (2.2), the 
points Pj influence all points in {M k (Afi))o.. UJ for sufficiently large k and for 
i — 1,2, respectively. Hence, < s x > s 2 for Sj = p ( j ) 1 (M^(A/'j)) and since 
Si influences all points in (M^Mi)o..j for all sufficiently large(r) k, we have 
derived Statement (b). □ 

(4.6) Lemma. (MIN-MAX relation) 

Let M {k) := (M*M) o.. w be a subdivided grid mesh with segment angle in 
(0, 7r). Then there exists a positive v such that for all k > 



MlN{M {k) )/MAX{M {k) ) > v . 

Proof. Let p be a maximum vertex in H(M {k) ), i.e., ||p|| = MAX(M {k) ). 

First, we apply Lemma (2.2) to obtain an r G N such that p influences 
every vertex q in H(M^ k+r ^). The refinement and averaging operators R 
and A preserve symmetry. Therefore, any vertex q in H(Ai^ k+r ^) can be 
computed by calculating midpoints of vertices or by projecting vertices or- 
thogonally into the spokes 5o and 60.5. According to Lemma (3.6), the first 
half segment lies in the cone C = C(M). Since C is non-negative, it follows 
for the midpoint of any two points a,b 6 C that 



(4.7) 



1 1. 
-a + -b 

2 2 



1 „ 

> - 2 M 
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Furthermore, if b is the projection of a into S0.5 and if c is the intersection 
point of So and the line ab, we have 

(4.8) ||b|| = W||c|| = cos(vV2) ||c|| > cos(p/2) ||a|| , 
||c|| 

as illustrated in Figure (4.9). 




(4.9) Figure. Projection into the spoke S0.5. 



Note that in deriving Equation (4.7), we used that the coordinate axes 
are perpendicular to the spokes S and S , 5 . Because of symmetry, the same 
holds if b is the projection of a into So. Altogether we see that there exists 
a constant /i > depending only on n and <p such that 

MIN(M (fc + r) ) > n MKX(M {k) ) . 

Second, we retrace the above argument with Equation (4.7) by 

max{||a||, ||b||} > I|| a || + ^||b|| 

1 1 
£ 2 a+ 2 b 

and Equation (4.8) by ||a|| > ||b||. This then shows 

MAX{M {k) )) > MAX(M {k+r) )) 
and thus 

Mm{M {h+r) ) > At MAX{M {k+r) ) . 
Because \i and r do not depend on k, Lemma (4.6) holds for 

v = min ({//} U {MIN(M^)/ MAX (-^ (fc) ) I k = 0, 1, . . . , r - 1}) . 

□ 
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5 Subdividing grid meshes 



In this section we show, that a sequence of subdivided grid meshes converges 
to an eigenvector of M n , where we use the results of Section 4 to prove that 
the sequence cannot have two distinct accumulation points. Let Af be a grid 
mesh with segment angle (p E (0, 7r) and j + 1 rings, where j > cu = [^y^J • 
Let Af {k) = (M*j\0o..j and Af k = M^ k) /\\M {k) || for k > 0. In this section, 
we show that the sequence of the normalized ringnets A4 k — A<o.. w /||A<o.. w || 
converges (Theorem (5.4)) to a symmetric positive eigennet of M n restricted 
to core meshes. 

Because of Lemma (4.4) (b), Mk is positive for all sufficiently large k and 
we can apply Lemma (4.6) showing that MIN(AAk) is greater than a positive 
constant for all sufficiently large k. Hence, and with the aid of Lemma (3.6), 
we obtain 

(5.1) Lemma. (/C > 0) 

Every accumulation point fC of (A4k) is positive and its first half segment lies 
in the cone spanned by the spokes S and S , 5 of M., 

H{K) c C(M) . 

Let M be the operator M n restricted to j-nets. It can be represented 
by a quadratic matrix with J\f^ = M k J\f because of Lemma (2.2). Let 
Ai, . . . , X p be the eigenvalues of M and let v°, . . . , be the eigenvectors 
and the generalized eigenvectors associated with Aj, i.e., 

Mv° = A iV ° and Mv[ = A,v| + v^ 1 , / = 1, . . . , on . 

With respect to these vectors, Af has a unique decomposition 

i=l 1=0 

with (7j < ccj and u it(Ti ^ 0. 

It can be readily checked that 

M k v\=^\ k - l v° + u k 

with ||ufc|| = o(| ( k )\ k ~ l \)- Assuming |Ai| = • • • = |A S | and (Ti = • • • = a s > 
and assuming for i > s that <jj < o~\ if |Aj| = |Ai|, or = —1 if |Aj| > |Ai|, 
we get the following lemma. 
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(5.2) Lemma. (Expression of accumulation points) 

The limit of any convergent subsequence o/(A4), i- e., any accumulation point 
of {Mk), lies in 

V:=span{v?,...,v°}. 

// Ai = • • • = A s > 0, then any accumulation point of (A4) is an eigenvector 
of M with eigenvalue X±. 

(5.3) Lemma. (Accumulation point 1) 

Here, we multiply and raise vectors to a power coordinate wise. The sequence 
{r k ) in C l with r = [r\ ... and \r\\ — • • • — \r\\ — 1 has the accumulation 
point 1 := [1 ... 1]*. 

Proof. The sequence (r fc ) is bounded and therefore has an accumulation point 
h. Hence, for all e > 0, there are integers l\ and l 2 with arbitrarily large 
differences l± — li > such that 

\\r h - h|| < e and ||r' 2 - h|| < e , 

where || • || := || ■ H^. Because |r' 2 | = 1, we have 

llr.il— '2 1 || Hfr.ll—12 ]_W 2 " 



lj.il j.h I 



< ||r Zl - h|| + ||r' 2 - h|| < 2e , 
which concludes the proof. □ 

(5.4) Theorem. (Convergence of (Mk)) 
The sequence {Mk) converges. 

Proof. Since {Mk) is bounded, it suffices to show that {Mk) has only one 
accumulation point. Let JC and C be two accumulation points of {Mk)- Using 
Lemma (5.2) with j = u, we can write /C and C as 

s s 

K, = z/jV° and C = /XjV° . 

i=\ i=l 

Let A = |Ai|, n = Xi/X and r fe = [rf ... r k }\ Then, 

y k M k JC = hv? . . . u.y s ]r k and ^M k C = ^v? . . . ^v°]r fe - 

Due to Lemma (5.3), there is a sequence {k n ) with r kn — > 1 as n — > oo. 
Therefore, 

(5.5) -^-M kn K -> K and v^M fe "£ -> £ . 
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From Lemma (5.1), it follows that /C > and C > 0. Therefore, there is 
a largest sequence (a k ) C M >0 and a smallest sequence C 1R >0 such that 



(5.6) a k 



M k JC < 



M k C < b k 



M K, 



for all k > 0. Since (a k ) is maximal and (b k ) is minimal, (ak) is monotonically 
increasing and (bk) monotonically decreasing because of Lemma (4.4) (a). 
Since (a k ) is bounded by b and (bk) by a , both sequences converge. Let 
a = limafe and b = limfo fe . Then we infer from (5.5) and (5.6) that 

aK < C < bK , 

and again using (5.6) that a < ao < bo < b. Consequently, the monotonicity 
implies = a and bk = b for all /c. 

Next, we assume C > aJC ^ C Due to Lemma (5.1), we can apply 
Lemma (4.4) (b) and get M k C > aM k K, for all sufficiently large k, which 
contradicts the maximality of = a. Thus, C = aJC. 

Since 1 = ||A^fc|| = = ||£|| = = a ll^l|) we obtain a = 1 and 

C = JC, which concludes the proof. □ 

(5.7) Corollary. (Properties of the limit ringnet) 

Let Aioo be the limit ringnet of (M.k)- 

(a) .Moo is an eigenvector of M with eigenvalue A = ||M«Moo|| > 0. 

(b) Moo > 0. 

(c) M.00 has norm 1. 

(d) Since M preserves symmetry, M.00 is symmetric with the same segment 
angle as M. 



Proof. Because 

MMoo 
\\MMoo\\ 



M lim fe 



\\M k M\\ 



M lim fc 



M k M 



M k+1 M 
fc™ ||M fe+1 A^| 



M k+1 M 
\\M k M\\ 



k— >oo 


M k +^M 




\\M k M\\ 



M c 



we get Statement (a), and by Lemma (5.1) Statement (b), and Statement (c) 
has been shown in the proof for Theorem (5.4). □ 
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6 Subdominant eigenvalues of the midpoint 
operator restricted to core meshes 

Let M be a grid mesh with u + 1 rings, frequency / and segment angle 
tp = 1f-r\fm G (0,7r). Let M be the matrix of M n restricted to core meshes 
and let = M k M for k > 0. Due to Theorem (5.4) and Corollary (5.7), 

the sequence Mk '■= M^ /\\M^\\ converges to a symmetric eigennet Moo 
with segment angle p and some certain positive eigenvalue A. In this section, 
we show that A is the dominant eigenvalue for the eigenspace of frequency 
/ and that A is smaller for higher frequencies. Moreover, we derive that for 
frequency 1, A is a double subdominant eigenvalue of M. 

(6.1) Theorem. (Maximum property of .Moo) 

Any rotation symmetric eigennet Af of M with segment angle <pj^ e [<p, ir] 
has an eigenvalue \i with < A or is a multiple of Moo- 

Proof. It suffices to consider the case that Af is not a multiple of .Moo- 

First, we assume that Af is reflection symmetric. As in Lemma (4.4), we 
rotate .Moo by (<pj^ — (p)/4: degrees counterclockwise. We use the coordinate 
system of Af and normalize Af such that 

Af < Moo and Af ft Moo ■ 

From Lemma (4.4) (b) it follows that for sufficiently large k 

M-"' A' < M 2k Moo and n 2k Af < \ 2h Moo ■ 

From Lemma (3.5) it follows that ji 2 is positive, which yields < A. 

Second, we assume that Af is not reflection symmetric and /i ^ 0. Using 
the scalar product (Af,Moo) '■= Af t M 00 , we define 



and obtain the symmetric ringnet 
Af s =Af 1 +AT 1 

with segment angle ip^f. Since {Af, Moo) = (Af, Moo), we have (Af s ,Moo) = 
and, hence, Af s ^ /3Moo forall /3 E C\{0}. 

If Af s = 0, then iAfi = iAfi is a symmetric ringnet with eigenvalue \i and, 
as shown above, < A. 




Af if (Af,Moo) = 
-Af if (Af,Moo)^0 
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If M s 7^ 0, we rotate Moo by (ipx — <p)/4: degrees counterclockwise as in 
Lemma (4.4). We use the coordinate system of M s and normalize M s such 
that 

(6.2) M 8 < and M s £ M^ ■ 

Because of Lemma (4.4) (b), there exists an a e (0, 1) such that for suffi- 
ciently large k 

/?M + Ji k M~i = M k M s < aM k Moo = a\ k M oc 

and 

\wJ \wJ \mJ 

According to Lemma (5.3), there is a sequence (/cz)/ £ n with {ji/\ji\) kl — > 1. 
Consequently < X because otherwise 

lim (-r^) M kl M s =M S < aMoo < Moo , 
i^oo y J 

which contradicts (6.2). □ 

The subdivision operator M can be represented by a block-cyclic matrix. 
Consequently, there is a basis of rotation symmetric (generalized) eigen- 
vectors with frequencies / e {0, . . . , m — 1}, see [PR98, Section 3], [ZS01, 
Section 4.2] and [Che05, Section 3.5 on Pages 47-58]. Since the conjugate 
eigenvectors have the frequencies m — f, we obtain 

(6.3) Corollary. (Dominant eigenvalues) 

The eigenvalue A of Moo is the dominant eigenvalue associated with frequen- 
cies f and m — f . (Since Moo an d A depend only on ip and the fixed n, we 
denote A by 

For m — 4, the eigenvalues X v are well-known. Around a regular vertex, 
the midpoint subdivision surface consists of four polynomial patches and 
subdividing these means to map any polynomial p(z) onto p(z/2). Hence the 
eigenfunctions of M are monomials and the eigenvalues are powers of 1/2. 
Since the least monomials with frequency and 2 are quadratic, besides the 
constant 1, we obtain that Hq — 1/4 is the subdominant eigenvalue for the 
eigenspace with frequency and that = 1/4. Furthermore, Theorem (6.1) 
implies 

(6.4) X a > Xp > X n = no = 1/4 for < a < (3 < n . 
Consequently A27r/ m is the subdominant eigenvalue of M. 
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(6.5) Theorem. (The subdominant eigenvalue of M n restricted to 
core meshes) 

The eigenvalue A 27r / m of M n restricted to core meshes is subdominant and 
has algebraic and geometric multiplicity 2. 

Proof. It remains to show that A = A 27r / m has algebraic multiplicity 2 or 
that there is no generalized eigenvector associated with A. We assume that 
T-L is a generalized eigenvector with frequency 1 and associated with A, i.e., 
MH = Xn + M^. Then 

H 8 = \(H + % 

is a symmetric generalized eigenvector of M with frequency 1 because 

MU S = [MU + M%/2 = (XH + Moo + XH + M^) /2 
= {XU + M 00 + xk + M 00 )/2 = XU s + M 00 . 

Therefore, it can be presupposed that T-L is symmetric. Because .Moo > 0, 
there exists an a > with aM.^ > %. The ringnet aA^oo ~ H is symmetric 
of frequency 1 and is non-negative. From Lemma (4.4) (a) it follows that 

(6.6) M k (aM 00 -H) > 

for all k > 0. Furthermore, we have 

M k (aM 00 -H) = aM k Moo - M k H 

= a\ k M oc - (\ k U + kX^Moo) 




< 



for all sufficiently large k, which contradicts (6.6). Therefore, the algebraic 
multiplicity of A is 2. □ 



7 Subdivision of larger ringnets 

As shown in Section 6, subdividing a core grid mesh with frequency 1 leads 
to the subdominant eigenvectors of M n restricted to core meshes. In this 
section, we generalize this result to arbitrary large grid meshes. 



19 



We consider a j-net M with j > to + 1, where the control points are 
arranged so that 

Mo. M 



M 



M b 

M c 

M u +2..j 



i e 



where M c consists of the corners of the ring A/L+i, i- e., M c = {pL+i u+i 
Z m } if M is primal or Af c = {pL+2,^+2 I * e ^ * s dual, and A/& consists 
of the other points in A/^+i- With this arrangement of the control points 
and according to Lemma (2.2), the matrix of M n restricted to j-nets has the 
lower triangular form 



S = 



A 

* B 

* * 

* * 



C 

* 



where 



(M n Af) .. w = A/V c 



0..w ) 







' A 




' Ao..oj " 


(M n M) b 




* B 







(MJV) . 



LU + l — 



A 

* 
* 



* 



C 



Ao. 



w+l 



Hence, the eigenvalues of S are zero or are the eigenvalues of the blocks A, 
B and C. 

(7.1) Lemma. (Spectral radii of B and C) 

The spectral radii ps and pc of B and C satisfy 

p B < 2~ n and p c < A~ n . 

Proof. Since C is non-negative, we get [HJ85, Corollary 6.1.5 on Page 346] 

PC < IICIU = IIC1IU, where 1 := [1 . . . I] 1 . 

The vector CI represents the corners of the (oj + l)-ring A/^ +1 , where Mo. M = 
0, Afu+2..j — 0, Mb = and M c = 1. One can easily verify that the vertices of 
these corners are < 4~ n , which concludes the proof of the second statement. 
The first statement can be proven similarly. □ 
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According to Theorem (6.5), Inequality (6.4) and Lemma (7.1), we get 



(7.2) Theorem. (The subdominant eigenvalue of M n restricted to 
j-nets for j > u + 1) 

The subdominant eigenvalue of M n restricted to j-nets with j > ui + 1 is 
^2n/m an d A 27r /m has algebraic and geometric multiplicity 2. 

(7.3) Theorem. (The characteristic mesh of M n ) 

Let M be a grid mesh with frequency 1 and with p + 1 rings, where p = 
|~|n — |], see Equation (2.3). Let S p be the matrix of M n restricted to p- 
nets. Let 

n k := \\S k M\\ and M k := S k M/n k . 

Then for n > 2, (Af k ) converges to a positive subdominant eigennet, called 
the characteristic mesh of M n . 

Proof. Let m k := ||(^jV) .. w || and let M k := {S k p J\f)o.. w /m k . According 
to Theorem (5.4) and Corollary (5.7) (b), the sequence (Mk) converges to 
a positive eigennet M^ with eigenvalue A = X2Tr/ m - The dominant and 
subdominant eigenvalues 1 and A of S p are positive, see Corollary (5.7) (a). 
Furthermore, since 

(7.4) n k >m k E Q(\ k ) , 

we obtain by Lemma (5.2) that (Afk) converges to an eigenvector Afc with 
eigenvalue 1 or A. If (Afoo)o..u were zero, then S p Afoo = fJ-Noo with \p\ < 2~ n < A 
according to Lemma (7.1) and Inequality (6.4) , which contradicts (7.4). Con- 
sequently (A/"oo)o..cj is a multiple of Moo and is positive and A is its eigenvalue. 
Using Lemma (4.4) (b), we have A fc A/"oo = M k J\foo > for all sufficiently large 
k and hence, A/"oo > 0. □ 

8 The characteristic map 

Let C be the characteristic mesh of valence m of the midpoint scheme M n . It 
defines the control mesh of a characteristic map, which is a spline surface ring 
consisting of m segments with 3 subsegments of \n/2\ 2 polynomial patches, 
as illustrated in Figure (8.2). 

(8.1) Theorem. (Characteristic map) 

The characteristic map of M n is regular and injective for n > 2. 
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(8.2) Figure. One segment and its 3 • 9 polynomial patches for n = 6,7. 



Proof. Let Q = [0,2] 2 \[0, l) 2 and let c(u,v) : f2 ->■ C be the first segment of 
the characteristic map rotated such that c(l, 1) > 0. According to [PR08, 
Theorem 5.25 on Page 107], the characteristic map is regular and injective if 
all partial derivatives c v {u,v) lie in 

Q = {x + iy\x,y > 0} . 

To show this, we consider a rotated grid mesh Ai such that the subdivided 
and normalized meshes Ai k = (M^A^)o..p/|| (M^Ai)o..p\\ converge to the char- 
acteristic mesh C. 

Let £ k and £ be the sets of all edge directions p° j +1 — P°j in the first 
segment of Ai k and C, respectively. These and other edge directions control 
the partial derivatives c v of the first segment. Furthermore, we add both Ui 
and iu to £ k and S, where u and ui are the edge directions of the spokes So 
and Si, respectively. Refining and averaging a mesh also means to average 
its edge directions by the masks shown in Figure (8.3). 




1/4 "~ ~ ^ 1/4 

(8.3) Figure. Masks for the edge directions. 



In particular, the directions in £ k are either, due to symmetry, parallel to 
ui and m , or obtained by iteratively averaging the directions in £ k -\ and 
multiplying these by positive numbers because of the normalization. Thus, 
we know that £ k lies in the cone spanned by £ k -i, i- e., the cone spanned by 
{ui,m }, and therefore, both £ k and £ lie in the larger set Q U {0}. 
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Moreover, since C0..1 is symmetric and has at most one zero control point, 
at least one of its edge directions is non-zero. Subdividing C, we can see that 
every element of £ is a linear combination of £ with non-negative weights 
and a positive weight for the non-zero element in the 1-ringnet. Hence £ has 
no zero elements and lies in Q. 

A partial derivative c„(u),u G [0,2] x [1,2] (or u e [1,2] x [0,2]) is a 
convex combination of directions in £ or in £ reflected at Si (or at iSq), where 
a reflected direction has a smaller weight than its unreflected counterpart. 
Hence, c v (Q) C Q as claimed. □ 

Theorems (7.2) and (8.1) show that Reif's C 1 -criterion [Rei95, Theorem 
3.6] is satisfied. So finally, we have established 

(8.4) Theorem. (C^-property of M n ) 

The midpoint scheme M n is a C l -subdivision algorithm for all degrees n > 2 
and for all valencies m > 3. 
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